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We report an experimental approach to entanglement characterization in high-dimensional quan- 
tum systems using Schmidt decomposition techniques. A particular experiment uses spatial degrees 
of freedom of biphotons. We present a technique to realize projective measurements in Schmidt 
basis, allowing us to measure the coefficients in Schmidt decomposition and to estimate the degree 
of entanglement directly. Issues of modeling the spatial part of biphoton amplitude with simple 
double-gaussian functions are discussed and shown to be in good agreement with experimental 
results. 
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I. INTRODUCTION 

Spatial entanglement in spontaneous parametric down- 
conversion (SPDC) was a subject of intense research dur- 
ing the last decade. Besides fundamental issues, spa- 
tial states of biphoton pairs offer a platform for high- 
dimensio-nal quantum states engineering motivating this 
interest. One can distinguish two complementary ap- 
proaches to spatial qudit engineering with biphotons: 
one using "pixel entanglement" and similar schemes 
0], and another one based on using high-order coher- 
ent (usually Laguerre-Gaus-sian) modes j5l-[la|. In both 
approaches achievable dimensionality and collection ef- 
ficiency are figures of merit. Dimensionality of effective 
Hilbert space is limited by degree of spatial entangle- 
ment. In pixel entanglement schemes, for example, the 
pixel size should be made smaller than the coherence ra- 
dius of the pump in the far zone, and since the pump is 
always divergent, even a plain wave, selected by point- like 
aperture would be correlated to a whole set of plain- wave 
modes in the conjugate beam. The same holds in general 
for other possible choices of modes. A natural question 
to ask is whether there exist a "preferred" basis among 
the multitude of possible coherent spatial modes, which 
consists only of pairwise correlated modes and allows us 
to fully exploit the spatial correlations in SPDC? 

The answer is, of course, positive and well known - the 
decomposition into a set of Schmidt modes has all the 
required properties. Since it was used for the first time 
by Law and Eberly (l7j . it has become a common tool 
for entanglement analysis of infinite dimensional systems 
in general, and of spatial states of photons in particu- 
lar. Nevertheless, direct experimental attempt to address 
spatial entanglement of SPDC biphotons in Schmidt ba- 
sis was not reported until the recent work of authors 
[l8| . Here we report a detailed description of experi- 
mental techniques used for Schmidt modes filtering and 
detection, models used to describe spatial entanglement 



in SPDC, and the structure of Schmidt decomposition 
revealed in our experiments. 

The paper is organized as follows: in Section II we 
briefly review the iTiain features of SPDC angular spec- 
trum with emphasis on spatial entanglement. Section III 
reminds the main ideas behind Schmidt decomposition 
analysis of spatial entanglement, explicit expressions for 
Schmidt modes in both Hermite- Gaussian and Laguerre- 
Gaussian basis are discussed in Section IV with emphasis 
on Hermite- Gaussian case. Section V gives a brief review 
of experimental approaches to mode transformations and 
some recent results on estimating the degree of spatial 
entanglement and, finally, the detailed description of our 
experiments with spatial Schmidt modes is given in Sec- 
tion VI. 



II. SPDC ANGULAR SPECTRUM AND 
SPATIAL CORRELATIONS OF PHOTONS 

Biphotons generated in the SPDC process have contin- 
uous frequency and angular spectrum. Let us consider its 
structure in some details. SPDC can be phenomenolog- 
ically described using the following effective interaction 
Hamiltonian (l9| : 

H= j (f'fx^^\r)El-\r)E'^+\r))E'^+\r) + Y{.c. (1) 



Here Ep is the classical amplitude of the pump field and 
E is the scattered field operator. Considering pump to 
be monochromatic, the first order of perturbation theory 
gives the following expression for the state of the scat- 
tered field: 



mc)+ / dkidk2^{ki,k2)\l)f,^ 
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where IS. ~ ki + k2 — kp, loi + U2 ^ ujp. In the case 
of coUincar phase-matching and under wide crystal ap- 
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Scanned detector 




of scanned detector's motion. One can also study two- 
dimensional distributions by changing the scanning di- 
rection. These distributions are related to the biphoton 
amplitude in the following way: 
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Figure 1: A simple experiment to observe correlations in 
SPDC angular spectrum. 



/ dki±.^ \'^{ki±x,k2±x)\' 



(6) 



If the biphoton amplitude may be factorized in two 
terms depending on the wavcvectors of each photon sep- 
arately: vE'(fcr±, fc2j_) = *i(fci_L) X 'i'{k2±), i.e. if the 
biphoton state is separable, both distributions coincide. 
If this is not the case, and the state is entangled, the 
width of conditional distribution AfeJ'^' is less then the 
^{ki, fca) = £p{ki± + k2±)J^{ki± - k2±), (3) ^j^th of unconditional one Ak[''\ It was proposed to 



proximation one can obtain the following biphoton field 
amplitude [l^-Ill: 



where £p{ki±^ + k2i_) stands for angular spectrum of the 
pump, and J-{ki± — k2±) is a geometrical factor deter- 
mined by phase-matching conditions. 

Authors of [111, [22| give the following expression for J^: 



^'(fci, fe) = JV£p{kij_ + fc2±)sinc 



L{k 



i± 



, (4) 



with L being the crystal length and A/" - a normaliza- 
tion constant. This expression is valid only in the case of 
small pump divergence. A detailed analysis shows that 
in the case of highly divergent pump one should take into 
account the dependence of refractive index for extraor- 
dinary wave on the propagation direction. This question 
was addressed in [J, |23|, where it was shown that for 
highly divergent pump experimental results are better 
described by the following biphoton amplitude: 
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(5) 



where kr, and 

^ ak^ 



are taken at kp± = 0. 



Spatial correlations in SPDC angular spectrum mani- 
fest themselves in significant difference of single-particle 
and two-particle (conditional) distributions obtained in 
experimental scheme of Fig. [T] 

The setup consists of two point-like detectors posi- 
tioned in the far field and a coincidence circuit. One 
of the detectors is fixed, the other one may be scanned. 
The distribution of photocounts of the scanned detector 
corresponds to unconditional single-particle distribution 
dw™'^°"'^ / dkii^x, while the coincidence distribution in 
such a scheme corresponds to a conditional distribution 
dw'^°"-'^ /dkx_\_x, where ki_\_x is a transverse wavevector 
component of one of the photons along the direction 



use the ratio of these quantities Rk = Afc^^-'/Afc^'''' as a 
quantitative measure of spatial entanglement (which was 
later referred to as Fedorov ratio) [23 | . Experiments have 
shown, that the value of Rk for SPDC biphoton states 
may reach values as large as Rk ~ 100, corresponding to 
high degree of entanglement [l|, 0] ■ 



III. ANALYZING SPATIAL ENTANGLEMENT 
WITH SCHMIDT DECOMPOSITION 

The most developed approach to quantitative analysis 
of spatial (and frequency) entanglement of SPDC bipho- 
ton states is based on using coherent modes decomposi- 
tion. Biphoton spatial state space is "discrctized" by 
switching from continuous distributions in plane-wave 
basis of the previous section to discrete distributions in 
a chosen basis of spatial mode functions ^i(fcij^.2_L)- For 
an arbitrary choice of mode functions in the decomposi- 
tion of spatial state for each of the photons, the biphoton 
amplitude takes the following form: 

oo 

^{ki±,k2±) = C^yC.(^i±)e,(fc2±). (7) 

i,j=0 

It turns out, that by appropriate choice of the basis 
mode functions one can transform the expression ([7]) to 
a single-sum form 



^(^li, k2i_) = J2 \/%Mkl±)i'^{k2±) 



(8) 



which is called Schmidt decomposition. In this case 
the basis functions "tpiiki^^) should be eigenfunctions of 
single-photon density matrix pi,2(^i_L7 ^i_l), arid Xi are 
the corresponding eigenvalues. It means that the ap- 
propriate mode functions may be found by solving the 
following integral equation: 



Pi,2{k±, k'^)-tpi{k'j_)dk'j_ = AiVi(fcj 



(9) 
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It is straightforward to note that for a factorized state 
only a single coefficient Aq in the Schmidt decomposition 
is nonzero. For a highly entangled state, to the contrary, 
Xi only slowly decrease with i. Since the biphoton state 
is pure, it is natural to use single-photon fon Neuman 
entropy as a measure of entanglement [l^ : 



A. 



(10) 



Calculating entropy is a rather demanding task, requir- 
ing the knowledge of all eigenvalues in the Schmidt de- 
composition. It is much more convenient to use " average 
number of Schmidt modes" determined by Schmidt num- 
ber as an operational measure of entanglement [? ]: 

i^=^J-^. (11) 

Degree of spatial entanglement of SPDC biphotons de- 
scribed by wavefunction Q was analyzed by Law and 
Eberly in [l3|. The pump was assumed to be gaus- 
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Authors derived an ana- 



lytical expression for Schmidt number by approximating 
— k2i_) function of ([4]) by a gaussian function; 

2 



1 



(12) 



where b 



is the waist of the gaussian function 

modeling J-{ki± — k2±). In the following we will call this 
procedure " a double gaussian approximation" . The value 
of Kg is determined by a single parameter ba and may 
reach very high values for 6(t ^ 1 and ba ^ 1. Numerical 
calculations performed for 



J'ikii_ — k2i_) = sine b{kij_ 



(13) 



showed that real value of K is larger than Kg for all val- 
ues of ba. Authors also proposed a method to increase 
the degree of entanglement by spatial filtering. It turned 
out that the number of significantly nonzero terms in the 
Schmidt decomposition increases if one cuts off part of 
angular spectrum with small values of k±. Influence of 
spatial filtering on Schmidt number was further investi- 
gated in [2^ , where authors study the effect of irises and 
finite detection apertures. 



IV. SPATIAL SCHMIDT MODES STRUCTURE 

The choice of mode functions set in ([7]) is quite ar- 
bitrary and may be determined by convenience in an- 
alyzing a particular physical situation. If the situation 
corresponds to a beam propagating in free space, it is 
natural to chose the solutions of paraxial wave equation, 
i.e. Hermite-Gaussian or Laguerre-Gaussian modes, as a 
set of mode functions. Moreover it turns out, that these 
functions form the Schmidt decomposition for SPDC 
with the gaussian pump of moderate divergence. 



A. Laguerre-Gaussian modes and orbital angular 
momentum entanglement 

Let us introduce polar coordinates {k±,(f>} on the plane 
of transverse wave- vector components k±. Laguerre- 
Gaussian modes are described by functions of the fol- 
lowing form: 



LGp/(fcl) (xLI,'I 



(Afcj 



exp 



kl 



2(Afc_L)2 



exp ( + i { 



(14) 



Here ip'(x) are associated Laguerre polynomials, Afcj^ 
is angular width of the mode. Radial index p corre- 
sponds to number of zeros in radial direction and az- 
imuthal index I - to the phase shift on the closed loop 
around fc_L = point (topological charge of the beam). 
One can show, that Laguerre-Gaussian beams with non- 
zero I possess orbital angular momentum (QAM) [27l.[28j. 
This is also valid for single-photon beams: a photon in 
Laguerre-Gaussian mode has 0AM of Ih. 

In Laguerre-Gaussian basis decomposition ([7]) takes 
the following form: 

^{ki±,k;^) = J2Y1 ci,\%LGp,i,{k^^)LGp,ijk;^). 

h,pi h,P2 

_ _ (15) 
It was shown in [23, [SO] that 0AM is conserved in the 
SPDC process, i.e. in the subspace of fixed pi and p2 
the identity Iq = h + I2 holds, where Iq is azimuthal in- 
dex of the pump beam.jisj This conservation law gives 
rise to 0AM entanglement of photons in the same sense 
as transverse momentum conservation leads to spatial 
correlations considered in Sect. HIl In the subspace of 
fixed radial indexes decomposition in terms of Laguerre- 
Gaussian modes takes the form of Schmidt decomposi- 
tion: 



/P1P2 



^'(fcli, 

^C^('«7')LGp,z(A;r±)LGp,(j„_,)(^2l). 



(16) 



Properties of such decomposition were studied in (33| 
where authors propose to use its width (which they called 
"quantum spiral bandwidth") to quantify the degree of 
0AM entanglement. Quantum spiral bandwidth esti- 
mates the effective dimensionality of 0AM Hilbert space, 
but it is defined only for subspaces of fixed radial indexes. 
To get rid of double sum and transform ()16|) to a Schmidt 
decomposition form one should use double-gaussian ap- 
proximation. 



B. Hermite-Gaussian modes as approximate 
Schmidt modes for SPDC biphoton field 

Another convenient set of modes which is somewhat 
less used in literature is the set of Hermite-Gaussian 
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modes: 



H„ 



exp 



fc2 



(17) 



2(Afc_L)2 



where H„(x) are Hermite polynomials, and {fca:,fcy} arc 
transverse wave-vector components. This decomposition 
was studied in [s^l- Authors investigate dependence of 
coefficients of the decomposition on modal structure of 
the pump (which is also expressed in terms of Hermite- 
Gaussian modes) and on the relative width of the gaus- 
sian function in (jl7p and the pump mode. In this case 
the spatial parity of the beam is the conserved quantity: 
parity of product ^Gnimi{kix,kiy) x HGnam^ (^'22:, ^2^) 
equals to the parity of pump mode. 

As well as for Laguerre-Gaussian modes, one can get 
rid of two indexes in decomposition (jl7p and transform 
it to a form of Schmidt decomposition using double- 
gaussian approximation HH] • Let us consider the case 
of moderately divergent pump, when one may not con- 
sider the linear term in ([5]) . |43| For small a we can make 
a substitution suyc{'^) exp(— 7^), where 7 is a co- 
efficient chosen to make both functions "close" to each 
other. A good approximation is provided by choosing a 
value of 7 = 0.86.(53] The biphoton wavefunction now 
takes the following form: 



1'(fci,fc2) oc exp(--^-ti^^|^2i^)cxp 



(fci± - k2i_f 
262 



(18) 

where a determines the angular bandwidth of the pump, 
and b — ^ Akp/^L - the phase-matching bandwidth. 
Since the wavefunction is a product of functions depend- 
ing only on ki^2x and only on ki 2y, it is sufficient to 
consider the problem in one dimension: 

'^{kix,k2x) = 
exp( 



2 , {,ki^ + k2xf. 



nab 



2a2 



I exp 



jkix - k2xf 
262 



(19) 



One can show, that solutions of © for such a wave- 
function have the form |3£ 



1pn{kij:^2x) 




(20) 



where (j)n{x) = (2"n!-y/7r) ^/2g ^^/2H„(a;). For corre- 
sponding eigenvalues and Schmidt number we obtain: 



A„ = Aab 



{a-bf 



(a -f 5)2(»+i) ' 



2ab 



(21) 



So we have the following form of Schmidt decomposition 
for SPDC biphoton state under the double-gaussian ap- 
proximation: 

*(fci,fc2) = 

^\/%A^'lpniklx)'tp7nikly) X Tjjnik2x)'4'mik2y). ^ ' 
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Figure 2: Phase holograms transforming a Gaussian beam 
into Hermite-Gaussian HGnm (a) and Laguerre-Gaussian 
LGpi modes (b). Color range corresponds to phase change 
from zero to tt. 



Degree of entanglement for this two-dimensional wave- 
packet is given by Schmidt number: 



K = X K, 



b'f/Aa^b' 



(23) 



Let us note once again, that (|22p is not the only pos- 
sible form of Schmidt decomposition. It may as well 
be described in terms of Laguerre-Gaussian modes as in 
[13, HI]. The value of Schmidt number is, of course, basis 
independent, as was explicitly shown in a recent preprint 
by Miatto et al. 
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In fact, the difference between 
these two representations corresponds to the choice of 
polar or cartesian coordinates on the plane of transverse 
momentum components. 



V. EXPERIMENTAL METHODS OF SPATIAL 
MODES TRANSFORMATION 

There are several well developed methods to transform 
a gaussian laser beam to higher Laguerre-Gaussian and 
Hermite-Gaussian modes. The general idea is based on 
introducing a spatially varying phase shift corresponding 
to the desired mode to a beam [38l - l40j . This is achieved 
when a beam is diffracted on a phase mask or a holo- 
gram, acquiring a desired phase structure. Examples of 
phase holograms corresponding to Hermite-Gaussian and 
Laguerre-Gaussian modes of lower orders are shown in 
Fig. [2] Usually these masks are combined with blazed 
diffraction grating, "highlighting" the first diffraction or- 
der, which allows to get rid of nondiffr acted components 
increasing the quality of mode transformation. Using 
this type of holograms allows effective transformation to 
relatively high order modes [4l| . 

In modern experiments phase holograms are usually re- 
alized using active liquid crystal spatial light modulators 
(SLM). Devices based on parallel-aligned nematic liquid 
crystals (VAN) matrices on silicon substrate (LCoS) are 
most commonly used for mode transformation tasks (4^ . 
Every pixel of such matrix works as a bircfringent phase 
plate introducing a phase delay for a polarization compo- 
nent parallel to the crystals alignment plane. The delay 
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Figure 3: Experimental scheme realizing projective measure- 
ments in spatial modes basis. PM - phase hologram, SMF - 
single mode optical fiber, L - lens coupling gaussian mode of 
the fiber with gaussian mode of the entrance beam. 

is determined by the orientation of the crystals and is 
controlled by voltage applied to the pixel. Commercially 
available devices have pixels of 10 fj.m size. 

Adapting the technique of mode transformations to 
the single photon level allows to realize both transfor- 
mations and measurements in spatial modes basis. A 
principal sketch of experimental setup for projective mea- 
surements is shown in Fig. [3] The crucial element here is 
a single mode fiber guiding only the fundamental gaus- 
sian HGoo = LGoo mode. This mode is coupled with a 
lens to a gaussian part of incident light, thus realizing 
a projection on gaussian mode. To realize higher order 
projections the hologram transforming a desired mode 
to a gaussian one should be used. In such a setup only 
the mode corresponding to a particular hologram passes 
through the fiber and gives a click in the single photon 
counter, thus realizing a projective measurement. 

As far as we know, this scheme was realized for the first 
time at the single photon level in the work of Zeilinger's 
group 0. Authors have experimentally demonstrated 
0AM entanglement of SPDC photons. Following works 
of the same group showed the possibility to engineer 
qutrit Q and higher dimensional qudit [43| states and 
to realize quantum key distribution protocols with spatial 
encoding. 

Experimental preparation of photons with 0AM and 
study of OAM-entanglement of photon pairs is a rather 
developed field with numerous experimental contribu- 
tions [a, 0, ITol - [l6j (see review [31 for a much more com- 
prehensive list of references). In context of our work it 
is irnportant to mention recent works of Boyd's group 
l45j . where the quantum spiral bandwidth was di- 
rectly measured. The same quantity was estimated with 
excellent precision using indirect measurements in (l5l |. 
Authors obtained the coefficients of decomposition in 
Laguerrc-Gaussian basis by Fourier transforming the de- 
pendence of Hong-Ou-Mandel-type interference visibility 
on azimuthal rotation of one of the interfering beams. 
Quantum tomography of states with 0AM was carried 
out in [4^, states up to / < 30 were realized, but 
only a subspace spanned by \—l, —I) , |— /) , \l, —I) , \l, I) 
was used for tomographical reconstruction. Same au- 
thors succeeded in preparing arbitrary states in Z = ±2 
subspace using superposition of phase holograms for 
Laguerrc-Gaussian and Hermite- Gaussian modes [l3l |. 

Up to our knowledge, no experimental attempts to an- 
alyze spatial entanglement in SPDC in Hcrmitc-Gaussian 



basis were reported. The question of properly adjusting 
the experimental conditions to make the detected set of 
modes as close to Schmidt decomposition as possible was 
not addressed as well. These were the main concerns of 
our work. 



VI. EXPERIMENTAL REALIZATION OF 
MEASUREMENTS IN SCHMIDT BASIS 

Our main goal was to demonstrate the possibility of 
performing direct measurements in spatial Schmidt basis 
and to experimentally demonstrate all features specific 
to Schmidt decomposition. We have chosen to work in 
Hermite-Gaussian basis. In this basis Schmidt decom- 
position is symmetric in indexes m and n, making it 
convenient for characterizing complete two-dimensional 
entanglement structure. 

First of all let us address the issue of applicability of 
double-gaussian approximation in our experimental con- 
ditions. We used a 2 mm BBO crystal pumped by a 
CW He-Cd laser with Xp = 325 nm wavelength. The 
crystal was cut for coUinear frequency-degenerate Type-I 
phase-matching. Angular bandwidth of phase-matching 
in such crystal (neglecting the pump divergence) - b pa- 
rameter in (|18p is 6 = 0.037. It was convenient for our 
purposes to select the value of pump divergence corre- 
sponding to a moderate Schmidt number. We focused 
the pump inside the crystal with a 150 mm quartz lens 
and measured the divergence - a parameter in ()18p to 
be a = (5.8 ± 0.1) x 10~'^.[5lj We calculated eignevalucs 
and eigenfunctions for the reduced single-photon density 
matrix, corresponding to the precise SPDC wavefunc- 
tion ([5]) numerically. The calculation was performed as 
follows: Hermite-Gaussian modes corresponding to the 
approximate function (|19p were chosen as a basis, we 
have restricted ourselves to 10 lower order modes (giv- 
ing the Schmidt number with 3 decimal digits precision) 
and calculated the matrix elements of the precise den- 
sity matrix in this basis. Diagonalizing the calculated 
matrix, wc obtained eigenvalues and eigenfunctions. The 
results for lower order modes are shown in Fig. 21 One 
can see reasonable correspondence, at least, we should 
expect that phase holograms for Schmidt modes should 
be close to those of Hermite-Gaussian modes of appro- 
priate divergence. We should note that measured waist 
size [53 of the pump beam in the focal plane of the lens 
was Wp = (25 ± 1) /im, corresponding to = 1.4. 
That means, the pump beam is aberrated and is not re- 
ally gaussian, and that may cause some deviations from 
Hermite-Gaussian shape of Schmidt modes as well. 

We have used an LCoS SLM with VAN matrix pro- 
duced by Cambridge Correlators. The matrix has 1027 x 
768 pixels of 10 fj,m size each. It is an 8 bit device, ca- 
pable of introducing a phase shift of up to 0.87r. Since 
larger phase shifts are required for our holograms we used 
a double refiection scheme. We used two polymer film 
polarizers in front and after the SLM to reduce the un- 
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Figure 4: Numerically calculated Schmidt modes (blue) and 
approximate Hermite-Gaussian modes of double- gaussian ap- 
proximation (magenta) and corresponding eigenvalues. 
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Figure 5: Experimental setup for transforming an attenuated 
laser beam to higher Hermite-Gaussian modes. SMF - single 
mode fiber; 01,2,3 - microscope objectives; L - 145 mm lens; 
P - polarizers; D - single photon counter. A phase modula- 
tor is shown as a transmitting mask, while in reality it is a 
combination of SLM and a dielectric mirror M (see text for 
details). 



wanted polarization rotations by an additional dielectric 
mirror necessary in such scheme (see insets in Fig. I5l9p . 

To check the quahty of mode transformation with this 
device we used the setup sketched in Fig. [51 We used an 
attenuated 650 nm diode laser as a source. The beam was 
mode filtered with single mode fiber and focused with a 
20 X microscope objective to obtain the divergence simi- 
lar to that of an HGqq Schmidt mode in Fig. |H and the 
waist at the position of the crystal. So we obtained a sin- 
gle mode gaussian beam modeling the zero order Schmidt 
mode of SPDC beam. The beam was coUimatcd with 
145 mm lens and after reflection from SLM was focused 
with a 8x microscope objective to a single mode flber 
followed by a single photon counter (Perkin Elmer) . The 
beam waist exactly coincided with the mode size of the 
flber ( 4/zm). We should stress, that we paid special at- 
tention to mode matching and optics were chosen in such 
a way that detection mode exactly coincides with calcu- 
lated HGqo Schmidt mode. Parameters of phase holo- 
grams were adjusted to minimize the detector counting 
rate, i.e. to ensure orthogonality of transformed modes 
to a fundamental gaussian one. We have actually ad- 
justed three parameters: the position of phase step for 
^Gio(oi) modes, which is determined by the beam po- 
sition at the SLM (in horizontal and vertical directions, 
respectively), and the distance between phase steps for 
HG20 modes determined by beam size at the SLM plane. 
These parameters deflne the shape of holograms for other 
modes in a unique way. If we define " visibility" for mode 
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Figure 6: Detector counting rate for an attenuated laser beam 
transformed to HGnm modes. The fiber tip is set to position, 
corresponding to 2: = in Fig. [71 



transformations as the ratio of counting rates with holo- 
grams for HGnm modes to that for untransformed gaus- 
sian mode: V = (i?oo - Rmn)/{Roo + Rmn), than for 
almost all of the modes with < m, n < 4 it exceeds 
97%, corresponding to high quality of mode transforma- 
tions. The visibility is slightly lower for HGqi and HGn 
modes (still exceeding 90%), which we explain by poorer 
adjustment in vertical direction due to the design of op- 
tomechanical components used. Histogram of counting 
rates for various modes is shown in Fig. [51 

To check whether the spatial structure of transformed 
modes is really close to Hermite-Gaussian, we scanned 
the fiber tip in the focal plane of 03 objective. The 
counting rate dependence on fiber position is determined 
by the convolution of a corresponding Hermite-Gaussian 
function and a fundamental gaussian mode of the fiber: 



R{x) (X 

J H„„(y2£/w)exp ^--^^ exp 



dx 



(24) 



where w is the gaussian mode waist. Experimental de- 
pendencies are shown in Fig. [3 and have a character- 
istic shape of double-peak curves. Distance between 
maxima depends on the mode number, and is shown 
in Fig. [51 for "horizontal" HGno and "vertical" HGa^ 
modes, together with theoretical predictions for Hermite- 
Gaussian modes. To plot the theoretical predictions cor- 
rectly we estimated the waist size by fitting the convo- 
lution for HGqo mode with a gaussian curve, obtaining 
w = (3.87 ± 0.07) /im. 

When the attenuated laser beam is substituted with 
SPDC radiation, the described scheme realizes projective 
measurements in Hermite-Gaussian basis. Full scheme 
of experimental setup is shown in Fig. [9l Pump was fo- 
cused to a 2 mm BBO crystal with a 150 mm quartz 
lens LI, a second lens L2 with F = 145 mm was set 
confocal with LI to collimate the beam. Pump radi- 
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Figure 7: Counting rate dependence on the position of fiber 
tip for an attenuated laser beam transformed to various 
modes. Fiber is scanned in liorizontal direction (1) and in 
vertical one (2). 
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Figure 8: Dependence of maxima positions for fiber tip scan 
on mode number for laser beam transformations. "Horizon- 
tal" HGnO (blue bars), "vertical" HGom (red bars) modes and 
theoretical predictions (solid black line and dots). Theoreti- 
cal predictions are calculated for gaussian waist a = 3.9 fim 
(see text for details). 



ation was cut off with a UV-mirror UVM, and SPDC 
radiation was frequency filtered with an interference fil- 
ter IF. We used filters with central wavelength of 650 
nm and bandwidth of 40 nm and 10 nm and did not 
observe any significant improvement of visibility for nar- 
rower filter. All the following results were obtained with 
a wide 40 nm filter. Photon pairs were split with a 50/50 
non-polarizing beam-splitter. An SLM was placed in the 
transmitted channel and after refiection the radiation was 
focused into single mode fiber placed in the focal plane 
of 8x microscope objective. In the reflected channel the 
beam was focused into the same single mode fiber or, 
alternatively, in a multi-mode fiber with core diameter 
of 50 fim. Multi-mode fiber allowed all spatial modes to 
be detected simultaneously. To transform modes in the 
reflected channel we used phase masks made of thin mi- 
croscope cover glass plates. This construction was used 
instead of the second SLM, although it has limited per- 
formance for broadband SPDC radiation and is only suit- 
able for modes with n,m < 1. Signals in both channels 
were detected by single photon counters connected to a 
coincidence circuit. 

If the doublc-gaussian approximation is valid and 



Figure 9: Experimental setup. LI - 150 mm quartz lens; L2 

- 145 mm lens; BBO - 2 mm BBO crystal placed in the joint 
focus of LI and L2; UVM - UV mirror cutting off the pump; 
IF - interference filter for 650 nm with 40 nm bandwidth; BS 

- non-polarizing 50/50 beam-splitter; 01,2 - 8x microscope 
objectives; PM - spatial light modulator (is shown as trans- 
mitting mask for simplicity, real alignment is shown on the 
inset); PM2 - phase mask made of thin glass plates; SMF - 
single mode fiber; SMF/MMF - single or multi-mode fiber 
depending on the experiment (see text for details); Dl,2 - 
single photon counters (Perkin Elmer) . A 200 vertical slit 
S was used in "ghost" imaging experiments. 



Schmidt modes are really close to Hermite- Gaussian 
ones, we expect coincidences to appear only when simi- 
lar modes are selected in both channels of the setup. If 
different modes are selected in the transmitted and the 
reflected channels, no coincidences should appear, since 
only terms with equal indexes are present in This 
fact may be used as an experimental criteria of how well 
Hermite- Gaussian modes detected in our setup approxi- 
mate precise spatial Schmidt modes. Fig. [10] (1) shows a 
histogram of coincidence counting rate for the case when 
HGmn with < TO, ri < 4 modes were consequently se- 
lected in the transmitted channel and HGqo mode was 
selected in the reflected channel. We obtained values of 
visibility V = (i?oo - Rinn)/{Rm + -Rmn), where Rynn is 
a counting rate for HGmn mode, over 94% for all modes 
except HGoi and HGn modes. Comparing the results 
to those of Fig. [BJ we conclude that non-100% visibility 
for these modes is rather a result of technical imperfec- 
tions, common for schemes with attenuated laser beam 
and SPDC, than some physical discrepancy between de- 
tected modes and Schmidt modes. So we can state, that 
our scheme indeed realizes projective measurements in 
Schmidt basis. 

As is expected, when a phase mask corresponding to 
HGio mode is inserted in the second channel, maximal 
coincidence rate appears for the same mode in the trans- 
mitted channel, which is clearly seen as a peak shift on 
the histogram of Fig. [10] (2). We believe lower quality of 
phase mask used in the reflected channel to be the origin 
of somewhat lower visibility in this case. 

Further evidence of similarity of Schmidt modes to 
Hermite- Gaussian ones may be obtained by analyzing 
the dependencies of single counts and coincidences on the 
fiber tip position in the focal plane of the focusing mi- 
croscope objective. We expect the dependence for coinci- 
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Figure 10: Coincidence counting rate for SPDC radiation for 
various masks in the channel with SLM. (1) - HGoo mode is 
selected in the reflected channel, (2) - HGoi is selected in the 
reflected channel. 



(1) (2) 




fiber tip position, x (^m) fiber tip position, x ((im) 

Figure 11: Coincidence (1) and single counts (2) rate depen- 
dence on the fiber tip position in the focal plane of the micro- 
scope objective in the channel with SLM for different modes. 
Fiber is scanned in horizontal direction. 

dences to be described by ([M)) . Experimental curves for 
the case when fiber in the transmitted channel is scanned 
are shown in Fig. [TTJ Note that the double-peak struc- 
ture characteristic for Hcrmitc- Gaussian modes appears 
only in coincidences dependence, while single counts be- 
have monoto-nously, as is expected for spatially multi- 
mode radiation. The maximal value of single counting 
rate, however, decreases with increasing value of mode 
indexes. We will pay special attention to this dependence 
later. Distance between maxima behaves analogously to 
the case of attenuated laser beam, as shown in Fig. [T^l 

We obtained same dependencies of coincidence count- 
ing rate when the fiber tip was scanned in the reflected 
channel (see Fig. [T^ . In this case single counts, obvi- 
ously, do not depend on the mode selected in the conju- 
gate channel at all. This effect is a straightforward con- 
sequence of intermodal correlations in SPDC and may 
be thought of as a sort of " ghost interference" [13] . We 
should note that almost zero coincidence counting rate in 
the central position of the fiber is an interference effect 
demonstrating spatial coherence of detected modes. So 
this result may be considered as an experimental demon- 
stration of one of the main features of Schmidt modes - 
their spatial coherence. 

One can also observe the shape of Schmidt modes di- 
rectly. For this purpose we substituted the fiber in the 



20-1 
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Figure 12: Coincidence counting rate dependence on the fiber 
tip position in the reflected channel (without SLM) for differ- 
ent modes. Fiber tip is scanned in horizontal direction. 




12 3 4 

mode number, m 



Figure 13: Dependence of maxima positions for coincidence 
distributions of Fig. [TT]on the mode number (red bars). Cal- 
culated dependence for Hermite-Gaussian modes with w = 
(3.0 ± 0.1) nm corresponding to HGoo waist size (grey bars) 
is provided for comparison. 



reflected channel with a multi-mode one, which serves 
as a "bucket" detector, collecting all spatial modes. A 
200 fim vertical slit was inserted in front of the objec- 
tive in the reflected channel. Dependence of coincidence 
counting rate on the slit position forms a "ghost im- 
age" of the mode selected in the conjugate (transmitted) 
channel. Experimental results for the first three modes 
HGoo, HGio, HG20 are shown in Fig.^^ Solid curves are 
fit with Hermite-Gaussian functions l53l. 




slit position (mm) siit position (mm) slit position (mm) 



Figure 14: Coincidence counting rate dependence on the 
slit position with multi-mode fiber in the reflected channel. 
"Ghost images" of the modes selected in the transmitted 
channel are observed. Results for first three modes: HGoo 
- (a), HGio - (b), and HG20 - (c) are shown. 
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Figure 15: Distribution of eigenvalues in Schmidt basis. (1) 
Single counts rate of detector Dl in the channel with SLM de- 
pendence on the selected mode. (2) Calculated eigenvalues of 
single-photon density matrix in Hermite- Gaussian basis (i.e. 
in a double gaussian approximation). 



As we noted above, single counting rate in the channel 
with SLM for the central position of the single mode fiber 
decreases for higher order detected modes. Indeed s state 
of the single photon of a pair is described by a reduced 
density matrix, which has the following form in the basis 
of Schmidt modes ipniniki): 



(25) 



where index i = 1,2 numbers photons. Counting rate of 
detector Dl: Rimn ~ {ipmn\ P li'mn) ~ >^mn IS propor- 
tional to the eigenvalue - " weight" of the corresponding 
mode in the Schmidt decomposition. Analyzing single 
counts rate in the channel with SLM we can determine 
the eigenvalues in Schmidt decomposition. Histogram of 
single counting rate for detector Dl is shown in Fig. [15] 
(1). It should be compared to the calculated eigenval- 
ues of single-photon density matrix in Hcrmite-Gaussian 
basis shown on Fig. [15] (2). We can use fidelity: 



(26) 



as a quantitative measure of correspondence between the 
measured and calculated distributions. Experimental es- 
timate for eigenvalues Amn^^ = (i?lmn - Ro)/ J2mn Rlmn 

takes into account normalization, implying unit trace 
for the density matrix, and substitution of a constant 
noise level Rq, originated both from stray light and non- 
diffracted SPDC radiation. Experimental data give value 
of fidelity F = 0.92 ± 0.03. 

Using directly measured eigenvalues we may estimate 
the Schmidt number. Let us analyze the one-dimensional 
section of eigenvalues distribution shown in Fig. [16] for 
HGmo modes. This distribution is excellently approxi- 
mated with a gaussian one: Rimo ~ Ro + Cexp(— an) 
(correlation coefficient R^ = 0.993). Using expression 
([^T|) for the Schmidt number it follows that A„+i/A„ = 
{Kx — l)/{Kx + 1), and it is straightforward to obtain 
the following estimate for the Schmidt number: 
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Figure 16; Normalized distribution of eigenvalues in Schmidt 
basis. Experimental results (red) and exponentially decreas- 
ing calculated values under double-gaussian approximation 
(grey). Both distributions are normalized to unit Schmidt 
number (i.e. divided by ). 



Using experimental results for HGom modes we estimate 
the "horizontal" Schmidt number to be A'^. = 3.1 ± 0.9. 
Analogous procedure for "vertical" modes HGon gives 
Ky 2.7 ± 0.5. Both values are in agreement with K = 
2.97 calculated from (|2T|) for experimental values of a = 
3.5 X 10"^ and b = 0.020 (values are taken inside the 
crystal). 

Experimental dependencies of eigenvalues are in good 
qualitative correspondence both with exponential behav- 
ior for double-gaussian approximation and with the re- 
sults of numerical calculations. Small qualitative discrep- 
ancies are still above the level of statistical fiuctuations. 
They can be partially explained by systematic fluctua- 
tions of signal level on long time scales caused primarily 
by temperature fluctuations. Experiments have shown, 
that both coincidence and single counts rates are ex- 
tremely sensitive to small beam displacements. So work 
of an air-conditioning system stabilizing the room tem- 
perature cau-sed small but noticeable periodic change of 
signal with a period of 10 min. Error bars shown on 
Fig. [16] correspond to the amplitude of this fluctuations 
since the exposition time was comparable with their pe- 
riod. Complete thermal isolation of the setup should re- 
move this source of errors. Another significant source 
of errors may be the difference between detected angu- 
lar aperture and full angular bandwidth of SPDC. Sharp 
dependence of counting rates on the position of fiber tip 
relatively to the microscope objective shows maximum 
in the position slightly different from the focal plane. 
Moreover, the pump, assumed to be gaussian, was ac- 
tually aberrated to = 1.4, which may also slightly 
change the real distribution. Nevertheless, we believe our 
results to be clear evidence of possibility to realize pro- 
jective measurements in spatial Schmidt modes basis for 
SPDC biphotons, and the implemented scheme is useful 
for quantitative measurements and state reconstruction. 
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VII. CONCLUSION 

We have analyzed spatial entanglement in SPDC in 
terms of spatial Schmidt decomposition and shown that 
under reasonable assumptions, applicable to the par- 
ticular experiment, these modes are close to Hermite- 
Gaussian modes. Using this fact we have experimentally 
realized a scheme of projective measurements in Schmidt 
basis using an active spatial light modulator. Experimen- 
tal results prove high quality of gaussian beam transfor- 
mation to higher order Hermite- Gaussian modes. For a 
spatial multi-mode SPDC radiation such spatial filtering 
allowed us to realize projective measurements in Schmidt 
basis. We have experimentally verified similarity of spa- 
tial Schmidt modes for SPDC to Hermite- Gaussian ones. 
We have performed direct measurement of eigenvalues 
in Schmidt decomposition and experimentally estimated 
the Schmidt number, obtaining good agreement with pre- 
dictions of a simple doublc-gaussian model of SPDC spec- 
trum. 

Spatial modes are interesting objects for quantum state 
engineering. The Hilbcrt space is vast and its dimen- 
sion may be chosen at will by simply focusing the pump 



properly. One can both obtain single-mode biphotons 
and highly entangled EPR-like states as extreme cases. 
However, intermediate case, considered in our work also 
offers interesting possibilities. We foresee at least two 
lines of research in this area: quantum tomography with 
high-dimensional spatial states and quantum state engi- 
neering of non-trivial states. Schmidt modes measure- 
ment techniques developed here should be a useful tool 
for both tasks. 

Another interesting question is the behavior of 
Schmidt-like intermodal correlations in the classical 
limit. What would stand for the Schmidt number, and 
whether the modal structure may be described in terms of 
some simple decompositions like in the case of SPDC, are 
interesting issues to address both theoretically and exper- 
imentally. Our results for quasi-thermal light sources will 
be reported elsewhere. 
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10-02-00204 and 12-02-00288. S. S. Straupe and I .B. Bo- 
brov are grateful to the "Dynasty" foundation for finan- 
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